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Abstract 

The spectrum of partons in a QCD jet becomes independent of the 
primary energy in the low momentum limit. This follows within the 
perturbative QCD from the colour coherence in soft gluon branching. 
Remarkably, the hadrons follow such behaviour closely, suggesting the 
parton hadron duality picture to be appropriate also for the low mo- 
mentum particles. More generally, this scaling property holds for par- 
ticles of low transverse and arbitrary longitudinal momentum, which 
explains an old experimental observation ("fan invariance" ) . Further 
tests of the perturbatively based picture for soft particle production 
are proposed for three-jet events in e'^e~ annihilation and di-jet pro- 
duction events in 7^, 77 and pp collisions. They are based upon the 
difference in the intensity of the soft radiation from primary qq and 
gg antennae. 
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1 Introduction 



The study of multiparticle production in hard processes provides one with 
an important information about both the QCD partonic branching processes 
and the transition from the coloured partons to the colourless hadrons. The 
vast amount of experimental data collected in various hard collisions pro- 
cesses convincingly confirms that the inclusive characteristics of QCD jet 
systems can be successfully described within the analytical perturbative ap- 
proach to multiparticle production, see e.g. [00,0 . This approach is based 
on the so-called Modified Leading Logarithmic Approximation (MLLA) 
and on the concept of Local Parton Hadron Duality (LPHD) ||^ . 

In particular, the experiment clearly demonstrates that the bright colour 
coherence effects survive the hadronization stage and are distinctively visible 
in the data. This has been known for quite a long time for e"'"e~ annihilation 
processes (see e.g. 0]). Recently the very impressive data from HERA 0,||J^ 
and from the TEVATRON |]10|,|lll have become available. 

The LPHD allows one to relate the (sufficiently) inclusive hadronic ob- 
servables to the corresponding quantities computed for the cascading par- 
tonic system. One of the well known (but still striking) predictions of the 
perturbative scenario is the depletion of the soft particle production and the 
resulting approximately Gaussian shape of the inclusive distribution in the 
variable ^ = ln(P/£^)for particles with energy E in a jet of energy P (the 
so-called "hump-backed plateau") |^,12,13 . 



This behaviour is observed in all inclusive energy distributions; the charged 
particle spectrum for not too small momenta prove to be in surprisingly good 
agreement with the MLLA-LPHD predictions reflecting the QCD cascading 
picture of multiple hadropro duct ion. Moreover, the data collected in various 
hard scattering processes {e~^e~ , DIS, pp) clearly demonstrate a remarkable 
universality of particle spectra assuming the proper (MLLA-based) choice of 
the cascading evolution variables, equivalent to the e~^e~cms energy y/s. 

The most challenging is the soft momentum end of the particle spectra 
(p < 1 GeV) where the non-perturbative dynamics could invalidate the per- 
turbatively based expectations. An attempt to stretch the perturbative pre- 
dictions to the limit of their applicability (or better to say, beyond it) has been 
performed in previous papers of the present authors (see e.g. [|T4| , [T5| , |T^ , |r^ ) . 
In these papers it was shown that (after the proper modifications) the per- 
turbatively based formulae allow a sufficiently smooth transition into the 
soft momentum domain. The non-perturbative hadronization effects are en- 
coded in the transverse momentum cutoff Qq which can be motivated by the 
space-time picture of hadroproduction in QCD jets, see e.g. [0. Explicit pre- 
dictions for distributions of hadrons in the soft region need some additional 
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assumptions on mass effects. 

Scaling phenomena are expected to be of more general nature. Let us 
recall that the gluons of long wave length are emitted by the total colour 
current which is independent of the internal structure of the jet and is con- 
served when the partons split. Applying the LPHD hypothesis one then 
expects that the hadron spectrum at very low momenta ph should be nearly 



independent of the jet energy 



As discussed in [|T5], the low momentum data support the basic ideas 
of QCD coherence and LPHD. Quantitatively, the analysis was performed in 
terms of the invariant particle density E-^ for e+e~annihilation into hadrons 
at low momenta in quite a wide cms energy region (from ADONE to LEP-2). 
The spectra were found to be in a good agreement with the scaling behaviour 
and with analytical perturbative expectations which became sensitive to the 
running of the coupling as at small scales [n.l^. The new HI data[Q follow 
these predictions as well, thus confirming the universality of soft particle 
production. 

In this paper we study further the phenomenology of the soft particle 
production for both the charged and identified hadrons. The analytical cal- 
culations are compared with the quantity ^ which has some practical ad- 
vantages over the invariant density considered before. We derive a scaling 
prediction for particles at low transverse but arbitrary longitudinal momenta, 
which generalizes our previous results and which for very small p± explains 
the so-called "fan-invariance" experimentally observed long ago |T9| . 

An important test of this line of approach towards soft particle production 
is the sensitivity to the colour charge of the primary emitters: the density of 
low momentum particles roughly doubles when going from a primary qq to 
a gg antenna. It is interesting to note that this prediction applies to the soft 
particles already at present energies, whereas for the total multiplicity this 
doubling needs much higher energies. The sensitivity of the low momentum 
particle density to the primary colour charges can be studied in experimen- 
tally accessible processes, such as in e'^e" 3 jets [Jl5[. Here this process is 
discussed further and we derive predictions for the soft particle production in 
quark and gluon exchange reactions which are observed in photoproduction 
or DIS at HERA or in 77 collisions. 
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2 Perturbative predictions for particle pro- 
duction in the soft limit 



In refs. [0,0 an approximate solution of the MLLA evolution equation 
for the inclusive energy distribution of soft gluons, originating from a 
primary parton A, has been derived. This solution can be written as 

D'A{^,Y) = ^D^^{^,Y)Uexp{-G) , G = ^ J%Mdy (1) 



with 



^I{^,Y)dl = P^I^ ^ drln^—ln ^— 

^ A JO A r + A 



(2) 



Here we have used the logarithmic variables ^ = ln(l/x) = ln(Q/£'), Y = 
InlQ/Qo) and A = \n{Qo/A) with E the particle energy and Q the jet vir- 
tuality {Q = PQ for a jet of primary momentum P and half opening angle 
6); Ca is the respective colour factor, i.e., Cg = Nc and Cg = Cp', 7o de- 
notes the anomalous dimension of multiplicity and is related to the QCD 
running coupling by 7o = ANcas/2n or 7q = ln(p_L/A) with = ANc/b, 
b = (llA'^c — 2nf)/3; A is the QCD-scale and Nc and Uf are the numbers of 
colours and of flavours respectively, a = ^-^c + The shower evolution 
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is terminated by the transverse momentum cut-off p± > Qo- 

The first term in (H) corresponds to the emission of a single gluon and 
yields the leading contribution for E ^ Qq. It is proportional to the colour 
charge of the primary parton. Furthermore, this term does not depend on 
the cms energy, contrary to the higher order terms which provide the rise 
of the spectrum for large E with increasing y/s. The spectrum vanishes at 
E ^ Qo as 

In order to relate parton and hadron distributions at low momentum 
P ^ Qo, one has to make some additional assumptions on how to include 
mass effects. In [|l4|,|l5l it was required that the invariant density Eh-^j^ for 
hadrons approaches a constant limit when p/j — >^ as observed experimentally. 
Here it is worth to recall that at low ph 

dn Wi{s,EhVs) 
d'^Ph s 
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where Wi{s, Eh\^) are the standard e+e analogues of the DIS structure 
functions Wi^q"^, z/), see e.g. pO|. As well known, Wi{s, E^y/s) are related to 



the matrix elements of the current commutators and should be regular when 
For spectra which vanish as in (|^) a plausible relation between partons 



and hadrons can be achieved by relating their distributions as [0)13 



Eh — 1 = KhEp — = KhD\[^E, y) (5) 

dph dpp 



at the same energies Ep = E^ = ypl + Qq = E > Qq where = = 
\n{Q/E) and Kh is a normalization parameter. If hadrons from both hemi- 
spheres are added, Kh should be replaced by 2Kh. For large energies Eh 3> 
Qo © coincides with the usual LPHD relation dn/d^ = KhD{^,Y). 

In |]15| we have considered the invariant density Edn/d^p = dn/dyd'^pT 
in the limit of vanishing rapidity y and transverse momentum p± or, equiv- 
alently, for vanishing momentum \p\ = p, i.e., 

^dn 1 ^ ^dn 
Io= lim E— = -\imE— (6) 

where the factor 1/2 takes into account that both e"'"e~ hemispheres are added 
in the limit p ^ 0. 

Then, since Eh — Qo ~ Ph/i'^Qo) small ph, the invariant hadronic 
density 

Eh^^=KhD\{i,Y)l{A'Kpl) (7) 
approaches the finite, energy independent limit as in (H) 



s-kNcXQI ■ 



^0 

In the fixed limit /?^/A is replaced by 70 and Iq ~ ^IQl- With prescription 
(|^) one obtains a good description of the moments of the energy spectrum 
Y) by the MLLA formulae with Qo = 270 MeV in a wide energy range 
||14| , as well as of the particle spectrum in the soft region {Eh < 1 GeV) |ll5| . 

Certainly at the moment there is no unique recipe allowing one to perform 
a smooth transition between the perturbative and non-perturbative regimes. 
Also there are certain approximations involved in the analytical results for 
low energy particles. In refs. [00] an alternative prescription has been 



discussed: 

dn f p 



d\np \E 



Dl{iE,Y)\i^s. (9) 
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where D\{^e, ^)|ls is the so-called Limiting Spectrum which solves the 
MLLA evolution equations for A = in the full kinematic region, except 
very close to the boundary E,e = Y. At this phase space boundary the Lim- 
iting Spectrum approaches the finite limit D^(Y,y) = which equals 
numerically 1.069 (1.055) for ?t,j=3(5). 

Basing on the prescription (^) one arrives at the alternative formula for 
the invariant hadronic density which in the case of e^e~annihilation reads 

with = j^D^g- As it is easy to see this prescription also guarantees 
that the hadronic density reaches a constant limit at ph 0. We shall 
compare here the expectation based on Eq. (p!0| ) with the experimental data 
for charged pions {Kh = Kt^). The MLLA Limiting Spectrum taken at 
A = Qo = 150 MeV has been found to provide a good description of pion 
spectra at relativistic energies, see e.g. @J2^. Recall that the low momentum 
region in charged particle spectra is dominated by the pions. 

3 Phenomenology of soft particle production 

The energy dependence of the soft particle production rate was investigated 
previously for e^e~annihilation in terms of the invariant particle density 
Edn/d^p in the wide energy interval (from ADONE to LEP-2) ||T5|. Re- 
cently, new data on soft particle production became available from DIS at 
HERA In the following we compare various results for the quantity 
dn/d^p as a function of the momentum p of the registered particle. When 
presented in terms of this observable the data points are independent of the 
choice of the theoretical parameter Qo; in this case, however, we are dealing 
with a quantity which is not directly related to the Lorentz-invariant phase 
space. Using this quantity may facilitate the comparison between different 
experimental results. 

3.1 Energy dependence of soft particle production 

Fig. 1 shows the momentum dependence of dn/d^p for various cms energies 
^/s in e"'"e~annihilation in comparison with the theoretical predictions. In 
Fig. la the distributions of charged particles are presented in comparison 
with the approximate MLLA result (|l|) with prescription (|]) for the relation 

''We are grateful to Yu. L. Dokshitzer for recovering this result. 
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between partons and hadrons. The theoretical prediction describes well the 
approach to the scaling limit for low momenta and also the shape and en- 
ergy dependence of the spectra with the choice Qq = 270 MeV, as already 
discussed in |jl5|. The detailed behaviour at the very small momenta p < Qo 
depends on the approximation at the parton level involving large production 
angles > Qo/p) and also on the assumed relation between partons and 
hadrons in this region. 

This point is illustrated in Fig. lb, where we show the theoretical predic- 
tions with the alternative choice of the low momentum behaviour discussed 
in Sect. 2, based on the approximate analytical MLLA solution for Qq = A 
(Limiting Spectrum) and relation @; this form describes the pion spectra 
for ^/s > 10 GeV with the mass parameter Qq = 150 MeV which is in agree- 
ment with the earlier fits to the pion spectra ||6|,|22|. With this prescription 



the particle production rate in the soft limit p — is lower than in Fig. la, 
but it also approaches an energy independent value in this limit and that is 
the essential prediction of the QCD approach based on coherent soft gluon 
emission. 

Looking at the data in Fig. 1 more closely, there seems to be a cer- 
tain increase of the soft particle rate at LEP energies. This is illustrated in 
Fig. 2 a and b where the data below and above ^/s = 91 GeV respectively 
in e"'"e~annihilation show separately a very good approach to an energy in- 
dependent limit at small p, whereas there is a difference in absolute rates 
of 20-30% in the two energy regions. This is also visible in Fig. 3a where 
we show the cms energy dependence of three values of the momentum p. 
We also include here the recent data from HERA [§] which closely match 
the e"'"e~results. At the lowest momentum p = 0.2 GeV we see an almost 
energy independent production rate, but a small jump at the LEP-1 energy 
is visible. A rise at the higher energies could partially be caused by the weak 
decay products of heavy quark particles (perhaps also particles with strange 
quarks) which would add incoherently to the spectrum of particles produced 
directly from the soft gluons. This could happen in particular at the be- 
cause of the larger branching ratio into 6-quarks. However, the experiment 
p3| clearly shows that at low momenta {p < 0.5 GeV) the effect induced by 
the bb events is below 5%. It would be interesting to clarify quantitatively in 
more details which part of the observed rise could be explained by the overall 
effect of the weak decays. A certain small increase of the rate is actually also 
expected from the theory at the nonzero momenta (see the curves in Fig. 
3a). For larger momenta the rise of the production rate with energy is more 
pronounced and this trend is reproduced by the theoretical calculation for 
the full range of energy where the data are available (a/s > 3 GeV). 

An approximate energy independence of the rate for the low momentum 
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kaons and protons is well tested for ^/s > 10 GeV where data at some 
common low momentum value are available. A comparison in Figs. 3a and 
3b suggests that the spectra become practically energy independent when 
the particle momenta are of the order of their mass or smaller. 

It should be noted that the energy independence of the soft particle rate 
is by no means trivial. We recall, in particular, that the particle density 
dn/dy at central rapidity ?/ ~ rises with cms energy. This is illustrated in 
Fig. 4 for TASSO and OPAL data where rather precise results are available 
at p = 0.25 GeV. Whereas in this energy range (14 to 91 GeV) the plateau 
height rises by about 70% the corresponding rise for the soft particles at 
p = 0.25 GeV is only about 10%. This is consistent with the QCD picture of 
particle production: the central plateau rises asymptotically with the same 
exponent ~ a/Iu s as the full particle multiplicity and this rise comes mainly 
from the high transverse momentum particles. The soft particle production 
is suppressed because of the coherence in the emission from all colour sources. 



3.2 Universal soft limit for different particle species? 

In Fig. 5 we consider for different identified particles with mass m the rescaled 
spectra nv'dn/d^p which are expected to approach energy independent di- 
mensionless numbers in the soft limit p — » 0. First we see again, that in the 
pion spectra the violation of energy independence occurs at lower momenta 
than in the spectra of the heavier particles. Furthermore, these rescaled spec- 
tra of pions and kaons, separated by an order of magnitude around p ~ 0.8 
GeV seem to converge towards a common value in the soft limit. A similar 
convergence is also suggested for the protons although the lack of precise 
data at small momenta exclude a firm conclusion. We find an approximate 
universality of the quantity 

rjfi 

C = m3 — Uo^O.1-0.2 (11) 

for the particle species considered here. This corresponds to our previous 
observation that the invariant density behaves like Jq = Edn/ d^p\p^Q ~ 
in the soft limit. In our theoretical picture this quantity is given by 
Jo ~ 1/Qo (and not 1/A^, for example). To the extent that Qq is related to 
the particle mass such a universality could be anticipated. It would be inter- 
esting to have more experimental information on low momentum particles to 
test better the validity of such a universal behaviour; new HERA data with 
identified particles at such small momenta could be helpful as well. 
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4 Scaling properties of particle production at 
small p_i 

4.1 Colour coherence and "fan invariance" 

As discussed before the gluon spectrum in the soft hmit becomes independent 
of energy of the primary parton above about 1 GeV. This energy indepen- 
dence can be observed, for example, in the cms of e"'"e~annihilation or in the 
Breit frame of DIS. This result can be easily generalized. We make a Lorentz 
transformation along the initial parton direction, then in the new frame the 
initial partons have different energies but the soft particle production rate 
should remain unaltered as long as the primary energies are still sufficiently 
large (and the colour of the primary parton has not changed). Consequently, 
the quantity 

Io{y)= l^E^= ^^^E-^ (12) 



should be independent of rapidity y for such transformations and this 
behaviour should still hold approximately for not much larger than the 
hadron masses. Therefore, if the coherence argument is applied to this class 
of frames one obtains the prediction of scaling for the distribution 

<fn 



f{y,P±)- (13) 



dydp^ 

Equivalently, for ?/ ^ 77 = — lntg('i9/2) with production angle one expects 
An approximate scaling behaviour of this type - but for the distributions 



normalized to unity - was found by the TASSO collaboration [|T^] in the 
energy range from 14 to 34 GeV long ago and was called "fan invariance" .0 
In the following we will show explicitly how this scaling behaviour arises in 
the perturbative QCD calculations; furthermore we predict how the above- 
mentioned scaling behaviour is violated. This will lead to the suggestion of 
other variables in which scaling is restored in a wider kinematical region. 



^We thank P. Mattig for pointing out to us this scahng law which triggered the study 
in this section. 
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4.2 Scaling properties of the analytical results 



The double differential distribution dn/dS,dd in ^ = ln(l/x) and angle d is 
derived from the ^-distribution of particles D{^,Y) in a jet of half opening 
angle by (see also |T^) 



d 



didlnd dln^""^^'' ^'®=''- ^^^^ 

with Y = \n{PQ/Qo) for small angles. For larger angles an appropriate 
variable is F = ln(2Psin(0/2)/(5o) instead of Y For the application 

at small p± we take the approximate expression (|Tp for D\{^, Y) and obtain 
(with = Ed) 



d'^n 



(16) 



c?^ d In 

Within the DLA this distribution (the first factor in the first term) is given 

by 



d^d\npj_ Nc 



1 



In 2^ 



+ 13^ In 



'In^- 



In 



'In^' 



(17) 



r - ^ + A = In 



y + A = In ■ 



In 



A' A A 

For small the first term in the brackets from the single gluon bremsstrahlung 
dominates and is energy independent. This is the term which corresponds to 
the "flat plateau" of radiation in the low p± limit; it is also responsible for 
the "fan invariance" as follows from the exact relation 

d^n d'^n , , 

(19) 



dp^di!) d^pdp± 

with = ln(P/p) defined in terms of the momentum p. The second term 
gives the energy dependent correction which vanishes in the soft limit p± 
Qq] the cms energy ^/s enters only through the variable P± = p±/x with 
X = 2E/^/s. The MLLA correction involves the exponent G and its derivative 



b{Y + A) 



(20) 



so it again depends only on F + A and ^ (for soft particles one can assume 
^ S> A). Therefore our perturbative results fulfil the scaling law 



d^n 
d^ d In px 



F{^,P±) 



(21) 
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i.e. the dependence on the cms energy enters only through the variable ^ or 
X. In this way the violation of the "fan invariance" can be partly absorbed 
into the scaling variables ^ or a;. 

The analytical predictions have been derived for the soft gluon radiation 
(with energies smaller than, say, 0{1 GeV)), while the angles are measured 
with respect to the primary parton axis. There are two effects which are 
beyond the perturbative calculations. 

First, the direction of the primary parton is not directly known, but may 
be related to the jet axis or to the initial hadron or photon direction. Here the 
situation can be improved by measuring the energy-multiplicity correlation 
PI taking the momentum of each particle in the event in turn as jet axis with 
the particle energy fraction as a weight. The angular correlation of this type 
has been studied using Monte Carlo events with favourable results ||36|] . 

Secondly, we can only perform calculations at the parton level, so in 
comparison with experiment an assumption on the effect of hadronization 
has to be made. This will be described in the next subsection. 



4.3 Relation between parton and hadron distributions 

According to the LPHD concept, we assume that parton and hadron spectra 
are proportional to each other if momenta are large compared to Qq. This 
picture has so far only been applied to observables which are integrated over 
p±, while we are now interested in distributions differential in p±. Since per- 
turbative formulae explicitly contain the cut-off > Qo? they cannot indeed 
be directly compared with the data at small pj_, but they have to be modified 
in order to provide a smooth behaviour in this region. There is no unique 
prescription to be followed in order to build this bridge between parton and 
hadron level, and what we are proposing here is a simple phenomenologi- 
cal ansatz which is consistent with kinematical requirements and the scaling 
behaviour in the soft limit. 

In analogy to the relation (|) applied to the energy spectra in Section 
2, we take the transverse momentum cut-off Qo as an effective mass of the 
hadron with = p1 + Ql and we then replace 

Pl,p=pl,h + Ql Ep = Eh (22) 

in eqs. (|T6D and ( p^Tf ) on the r.h.s. as well as on the l.h.s. This choice guar- 
antees the proper kinematical limits pj_^p > Qo and p±^h > for parton and 
hadron spectra, while they coincide for large p±. Furthermore we use as 
before E^\h = -^\p- So for hadrons one can write (dropping the label h 
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everywhere) 



with the function F defined by (0) at the parton level, whereby the replace- 
ments (^) are applied. Once again, relations such as ( p2D are not an essential 
part of the perturbative LPHD picture, rather they should be considered as 
a plausible extension into the region p± < Qq, consistent with kinematics 
and scaling properties. 



Combining now (|23D and one obtains quantitative theoretical pre- 
dictions for the angular distribution of particles. Results for different values 
of p\\ at ^/s = 14 GeV are shown in Fig. ^ for Qq = 0.5 GeV and A = 0.05. 
Choosing such parameters, the results describe reasonably well the experi- 
mental data obtained by the TASSO collaboration in the region of low 
P|l , while some deviation is visible for larger values of pyQ. The large values of 
P|| and of angles correspond, however, to large values of the transverse mo- 
mentum p±, where our approximations are not valid any more. This point 
is particularly important for the study of "fan invariance", i.e., the energy 
independence of the renormalized angular distribution: the distributions pre- 
sented by the TASSO collaboration are normalized to unity in the full 
angular region and therefore not directly accessible by our calculations. 

A better test of the scaling law (|T^) can be obtained from the direct study 
of the double differential distribution d'^n/dp\\dp± in absolute normalization 
(i.e., normalized to the particle multiplicity) as a function of the transverse 
momentum p± at fixed p\\. Fig. ^ shows the theoretical predictions in MLLA 
with the same parameters as in Fig. ^ at p\\ = 0.8 GeV at four different 
values of cms energy ^/s. A good scaling behaviour is visible in the region of 
very low transverse momentum, whereas the position of the peak and the tail 
are found to smoothly depend on the cms energy. If the distributions were 
rescaled to the same area as in [1^, the overall energy dependence would be 



rather small. 

In order to better understand the origin of this energy dependence, it is 
useful to separately look at the contribution of the different terms in the the- 
oretical formula. This is shown in Fig. p, where the d'^n/dp\\dp± distributions 
at p\\ = 0.8 GeV and ^/s = 34 GeV are shown for different approximations: 
the DLA (Born term only), DLA (Born plus next order term) and MLLA. 

^This value for Qq is largely determined by the position of the maximum in the angular 
distribution and is larger than the value taken above for the energy spectrum. A unification 
is possible in principle by calculating the energy spectrum from the integration over the 
modified double differential spectrum cPn/dEdp±, but this is not attempted here. 
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Whereas the Bom term of DLA, which is independent of cms energy, is 
dominant at small pj_, the other terms become important with increasing p_i_ 
and give rise to the energy dependence visible in Fig. |^. As the next order 
term becomes important for larger p^, we also compare in Fig. ^ with the 
prediction from the Limiting Spectrum which corresponds to an all order 
summation within a high energy approximation (taking the same Qo, ^ = 
and normalization {Kh = 3.5) adjusted to match the MLLA result). It is 
remarkable how close are the predictions from the MLLA and from the Lim- 
iting Spectrum. This result suggests that retaining only the first two terms 
in the DLA solution gives indeed a reasonable approximation in the region 
of low transverse momentum. 

The deviation from the scaling behaviour (|T^ is due to the weak log- 
arithmic dependence of the correction terms on the primary energy. As a 
consequence of ( pSD with p ~ py the invariant density, after our modifica- 
tion (^), still follows the scaling behaviour 

E^ = fixE,P±) , (24) 

i.e., for fixed p± it depends on particle and cms energies only through xe = 
2E / ^/s. This scaling behaviour (^4]) is expected to hold in a wider kinematic 
region than eq. ([T^). 

5 Further tests of the perturbatively-based 
picture of soft particle production 

5.1 Sensitivity to the colour charge of primary partons 

The nearly energy independent soft end of the particle spectrum has been in- 
terpreted above as a consequence of the coherent soft gluon emission together 
with the LPHD. Alternatively, one may argue that at these low momenta of 
few hundred MeV we observe the purely hadronic phenomena not related to 
the perturbative QCD at all. It is therefore important to test further the 
perturbative nature of the soft particle production. 

The production rate in the limit (H) cannot be predicted in absolute 
terms, but we note the dependence of this rate on the colour charge factors 
Cp and Nc for the jets originating from a primary quark or gluon respectively. 
Therefore, we have argued in that the sensitivity of the soft end of the 
spectrum to the colour of the primarily produced hard partons provides a 
crucial test of the perturbatively based picture and can teach us about the 
region of its validity. 
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The prediction for the ratio of multiphcities in gluon and quark jets Rg/q = 
Nc/Cp = 2.25 has been made long ago but holds only at asymptotic 
energies. At existing energies this ratio is measured to be smaller but not in 
disagreement with the perturbative QCD analysis (see, for example [^,^,0). 



Our present discussion concerns only the soft particles: in this kinematic 
regime the asymptotic predictions on the spectrum, in particular the scaling 
properties, are well satisfied and therefore the asymptotic ratio Rg/g in this 
limit could be asymptotic too without contradicting the global results. 

Therefore, in order to test this scenario, one would like to compare the 
soft particle production in qq and gg events. Whereas there are well measured 
spectra down to small momenta of about 200 MeV from quark jets in e"'"e~and 
DIS processes (see Figs. 1 and 2) the direct production of a gg final state is 
more difficult .[] To perform such a comparison nevertheless we have discussed 
I3| two different approaches. 



m 



First, there is the comparatative study of two jet events arising from 
reactions with primary hadrons and photons mediated by quark or gluon 
exchange processes. If the exchange is sufficiently hard the two types of 
processes yield primary colour sources of the desired qq and gg types. We 
will dicuss below how these limits can be derived explicitely from the hard 
scatterings in the limit of small momentum transfer. The phenomenological 



studies along this line are described in papers []T3 . 

Secondly, one may consider final states with more than two partons. In 
certain coUinear limits such final states approach a qq or gg type antennae. 
A simple example is provided by the process e'^e~ — > qqg which for the 
configurations with angle Qgg ~ or Qgg ~ behaves like a qq and for 
©59 ~ behaves like a gg final state [^]. Such extreme configurations are 
in general not realistic. However, we can study the radiation perpendicular 
to the event plane (say, in the cone with not too large opening angle 6), 
then the particle density at low momenta should vary by the factor Nc/Cp 
when going from the qq to the gg type configuration and this variation can 
be predicted also for the intermediate kinematical range. Similarly, one can 
study such an associated perpendicular radiation in parton-parton scattering 
processes with quark or gluon exchange at any angle (see Fig. 9). 

In the following we will present the expectations for various processes; 



the main results for e^e — > qqg have been given already in 



^At high energy e+e cohiders the process 77 gg may become a prospective source 
of clear information Ell 
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5.2 Soft radiation in e+e qqg 

It was realized long ago that the overall structure of soft particle angular 
distribution in multijet events in hard processes is governed by the underlying 
dynamics of colour at small distances. An instructive known example is the 
perturbative explanation |^ of the particle angular flow in e~^e~ qqg 



events ("string effect" 0]). The multiplicity flows build up a colour portrait 
of an event which can be used as a natural "partonometer" mapping the 
primary interaction short-distance process PD| , |55| ,P|. In these studies the 



drag effect on the particles within the production plane is caused by the 
primary colour dipoles. More recently, the azimuthal angular distribution 
around the jet directions 0j51| have been studied. The partonometry ideas 
are strongly supported by a wealth of experimental data from e~^e~ — *■ qqg 
annihilation (see e.g. |^) and the TEVATRON pp collider [p!0| . 



In the present application we are interested in the particle production at 
low momenta p of order Qq. A straightforward calculation is possible for the 
radiation perpendicular to the production plane of the jets. Otherwise, in 
multijet events, the cut-off condition p± > Qo is more difficult to implement. 
Also, in this case p± appears as argument of and does not change for boosts 
of partons within the plane. In the coUinear limits for the 3-jet events one 
obtains the radiation transverse to the jet for which the coherence arguments 
apply in the first place. 

In the simplest case of a boosted qq pair, or equivalently, in e+e^ qq'j 
the radiation pattern of gluons with momentum p in direction n is given by 
PI 



dNgg a.. 



dQfidp (27r)^p ^ ^ 

W^^in) = 2CAfj), (0) = ^ (26) 

aiUj 

where a^j = (1 — Uiflj) and = (1 — nn^) for q and q in direction and fij. 
For the radiation perpendicular to the plane we simply obtain 

(ij) = 1 - cos Qij (27) 

with the relative angle Gjj between the primary partons. This yields the two 
limits (ij) =0,2 for parallel and antiparallel momenta respectively. 

As we are interested only in the soft particles and the relative variations 
between different processes or configurations we restrict ourselves to the sin- 
gle gluon bremsstrahlung of order ag as in (|25|). The MLLA corrections will 



become important for larger momenta (see e.g. the discussion in ||43| , [15| 
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The soft gluon radiation in the qqg process is given by the antenna pattern 
(defined as in ([25|) ) 

W^^^in) = Nc[{l+) + (M - ^(+^)] (28) 

where (+, — , 1) refer to (g, q, g) and this bremsstrahlung pattern holds, irre- 
spective of whether the partons are incoming or outgoing. With (^) one can 
obtain the particle density in the direction perpendicular to the production 
plane. 

It is convenient to normalize the production rate of the soft particles to 
the corresponding rate in the standard process e^e" — > qq in its rest frame 
with W'^ = 4Cf according to ( ]26|) and define the ratio 

- dNf/dn^dp ~ aCf' ^ ' 

for a general process p. For the perpendicular radiation in the 3 jet events 
one finds 



RT{V) = ^[2 - cos 91+ - cosei_ - ^(1 - cose+_)]. (30) 

In our approximation the momentum spectrum does not depend on the angles 
Qij between the jets. It will be interesting to study the spectrum dn/d'^p in 
such a cone experimentally and to find out whether the angular dependence 
of the particle density in ( PPD holds down to small momenta. 

It is also interesting to note the difference of this prediction to the large 
Nc approximation in which the qqg event is treated as a superposition of 
two qq dipoles (see, e.g., |0). In this case the last term in eq. (^) drops 
out and Cp = (N^ - l)/2Nc ~ Nc/2. This yields 

^[2-cos0i+-cos0i_] (large A/'c) (31) 

Predictions from these formulae for Rx are presented in Table |l| for various 
relative angles Qij. Note, in particular, the limiting case i?^ = 1 for soft or 
collinear primary gluon emission and the proper gg limit for the parallel qq 
(0+„ = 0) configuration, as expected. 

The role of the large- iVc limit can be investigated also by studying the 
production rate in 3-jet events normalized to the sum of rates from the corre- 
sponding 2-jet events (dipoles) with opening angle 9i+ and 9i_ respectively: 



_ dNV' 

Ni f 1 1 - cos(ei+ + ei_ 



" dNf {Q,+) + dN fiOi^' ^^^^ 



iV^-lV iV^2-cos0i+ -cos0i_ 
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R± 


(large Nc) 


0i+ = 7r-0i_ 


1 


1 


(collinear or soft gluons) 

01+ = 0+_ = fTT 


1.21 


1.18 


01+ = 0+_ = fTT 


1.42 


1.35 


01+ = 0i_ = fvr 


1.59 


1.5 


(Mercedes) 
01+ = 01- = vr 
{qq antiparallel to g) 


^ = 2.25 


2 



dNfydNf from 



OD and its large- 



Table 1: Prediction for the ratio i?j 
A^(7-approximation (31) for different configurations of the qqg events (0i+ 
0„„, 0i_ = 0„,-, 0+_ = 2n- 01+ - 0i_). 



This ratio measures directly the deviation from the large- A^^c limit R±= 1 
and thereby from the gg-dipole approximation. This approximation is not 
necessarily limited towards soft particle production. For the simple case of 
Mercedes-like events (0i+ = 0i- = 0+-) one obtains R± = 17/16 = 1.06. 
The 2-jet rates for relative angle Qij which appear in the denominator of 
eq. ( |32D could be found experimentally from the corresponding qq'y final 
states. 

Let us list a few further results: 

a) A particularly simple situation is met for Mercedes-type events where no 
jet identification is necessary for the above measurements. 

b) The large angle radiation is independent of the mass of the quark (for 
> mQ/Ej^t)- 

c) The above predictions for the ratios R±,R± are derived for the soft par- 
ticles according to the bremsstrahlung formula (^5]) . One may also con- 
sider the particle flow integrated over momentum, as in the discussion of the 
string/drag effect. In this case one has to include all higher order contri- 
butions which take into account the fact that the soft gluon is part of a jet 
generated from a primary parton. Then the angular flow dN/ dVLfi is given by 
the product of the radiation factor W{n) and a "cascading factor" [Q. For 
the ratio of multiplicity flows one obtains the same predictions (^) and (|32|) 
as the cascading factors cancel. It will be interesting to see to what extent 
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the predicted angular dependence for both quantities - the multiphcity flows 
and the soft particle yields - are satisfied experimentally. 
d) The similarity of particle flows in 3-jet and 2-jet events with corresponding 
angles 6i+, Oi-, as expected in the large-iV^ approximation, should also 
apply to further details of the final state such as the particle ratios. Since 
the Lorentz transformation along the boost direction produces a larger drag 
for heavier particles, one may expect that after the boost the K/ir and p/n 
ratios for soft particles decrease and then the same is true for soft particles 
in 3-jet events in comparison to 2-jet events in their rest frame^. 



5.3 Soft radiation accompagnying photo-production of 
di-jets 

5.3.1 Direct and resolved processes 

Next we consider the photoproduction of di-jets. The soft bremsstrahlung 
depends only on the momentum vectors of the partons but not on the vir- 
tuality of the photon. Our results can be applied to jp as well as to 77 
collisions. The di-jet photoproduction has been studied recently in detail at 
HERA |]46| , ^ . In di-jet production one can distinguish in the leading order 
QCD approach the direct and the resolved processes In the first case 



the photon participates directly in the hard scattering subprocess and trans- 
fers a large fraction [x^ ~ 1) of its primary energy to the secondary jets; in 
the second case, the hard scattering subprocess involves the partons in the 
protons and also in the photons and the energy fraction is smaller {x-f < 1). 
The data indeed show two peaks in the distribution of x^ corresponding to 
the two types of processes. 

The direct processes are mediated by 

a) l{Pi)+9{P2) q{P3) + q{P4) (33) 

b) i{pi) + q{p2) q{p-i) + g{Pi), (34) 

the so-called photon-gluon-fusion and QCD-Compton processes respectively. 
Their differential cross sections are given by 

da" TTaasQl , ^\ for^\ 

da^ SnaasQl ft s\ 

-df = ^;^[-s^-t) ^^'^ 

^This expectation has been verified in the JETSET Monte CarloU]. We thank T. 
Sjostrand for providing us with this information. 
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where Qq is the quark charge and 



s = {pi+ J92)^ t = {Pi- Psf and u = {pi - p^f (37) 
t = -s(l -cose,)/2 and m = -s(l + cose,)/2 (38) 

with the scattering angle in the di-jet cms. The cross section is maximal 
for small t (or u) where it is dominated by quark exchange. 

In the resolved processes the partons (g, q and g) in the photon interact 
with the partons of the proton. We restrict ourselves in this discussion to 
the small angle scattering where the gluon exchange contribution dominates 
and consider the subprocesses 







9iP3) 


'+9iPi) 


(39) 


d) 


g{pi) + q{p2) 


9{P3) 


+ QiPi) 


(40) 


e) 


q{pi) + q'{p2) 


(l{P3) 


+ q'iPi) 


(41) 


f) 


q{pi) + q{p2) 


qips) 


+ (i{Pi) 


(42) 


9) 


QiPi) + q{P2) 




+ q{Pi) 


(43) 



(in processes d and / the g's could be replaced by the g's). 

The different exchanges in the direct and resolved processes cause different 
distributions dn/dt ~ and dn/dt ~ respectively and this difference 
is clearly seen in the data Therefore in this case, by choosing suitable 
intervals of x-y, one can select processes of different colour structure and 
one can check whether or in which kinematical interval the low momentum 
particles are sensitive to it. To this end we calculate the respective production 
rates for the soft particles. 

For the given total energy the hard process depends on 3 variables which 
we choose here for simplicity as the incoming parton energy fractions Xi and 
X2 and the momentum transfer t. The cross section for producing the di-jets 
with and without accompagnying soft gluons is then given by 

MSS* - E (44) 

dcj dcT^^^^^ 

, , = E/*(^i'/^')/i(^2,/i')— 1— (45) 
dxidx2dt dt 

where fi{x,fi'^) denotes the parton structure functions at scale /x^. We are 
interested in the density of the soft radiation, i.e. in the ratio 

dN da , da , , . 

- (46) 



diljfdp dxidx2dtdVtfidp dxidx2dt 
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where often in practical applications an integration over ranges of the kine- 
matical variables Xi,X2 and t has to be performed. If only one subprocess 
[ij kl] contributes we come back to relations analogous to 



5.3.2 Soft particles in direct processes 



The angular factors required in (p| 
in processes a) and b) in (|33|) and 
relabeling of the partons (see also 



W(n) 



c 



for the partons q 
41) are given by 



(23) + (24) - ^(34) 



N 



c 



c 



(34) + (24) - ^(23) 



c 



q and g involved 
with the proper 



(47) 
(48) 



For the case of perpendicular radiation in the di-jet rest frame we have simply 

(23) = (14) = 1 + cos e„ (24) = (13) = 1 - cos e„ (12) = (34) = 2. 

(49) 

If either of the two direct processes a) or b) dominates we obtain for the 
respective rates, normalized as in (B^) 



R'. 



1, 



c 



AC, 



3 — cos O, 



^(1 + cose,) 



(50) 
(51) 



This means that the soft radiation in the photon gluon fusion process is the 
same as in the standard process e^e~ —>■ qq irrespectively of the scattering 
angle Qg. In case of the QCD-Compton process we find in the limit of small 
angles ^ 1 which is dominated by quark exchange = 1- In the case 
of backward scattering with 0^ ~ vr (i.e. the transition 7 — (yf at small 
angle) we have effectively a colour octet channel and find R^j_ = Nc/Cp as 
in case of a gg final state. In the large Nc approximation we can consider 
the bremsstrahlung in the quark exchange processes as coming from a single 
colour dipole whereby the second is inactive (for = 0); in case of octet 
exchange the radiation comes from two dipoles and adds incoherently. 

If the two jets are not identified one has to add the cross sections for 
scatterings at angles 6^ and tt — Gg. One finds 



R 



a,sym 



R 



b,sym 



Nc G{x) [1 + x- y/Nl] + G{y) [I + y - x/N^] 
2Cf G{x) + Giy) 



(52) 
(53) 
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with X = —t/s, y = —u/s and G{z) = z + z ^. So the latter ratio rises 
from i?^'^™ = 1 at e, = to i?^'^"" = H = i.63 at 9, = f . Therefore, in 
the kinematic region where one of the processes dominates one can check the 
variation of the soft particle density with scattering angle B<j, analogously 
to the case of e'^e" qqg. In the more general situation one has to use the 
superposition ( ^41) with the structure functions as weights. 

In the further discussion we restrict ourselves to the scattering at small 
angles where quark or gluon exchanges dominate to simplify the discussion. 
Then in both processes a) and b) the quark exchange with R^^ = 1 dominates 
and this result also holds in a general superposition ( PD of symmetrized 
direct processes. 



5.3.3 Soft production from resolved processes 

For small angles —t<^s the expressions for the radiation patterns simplify 
significantly I^H]. The antenna patterns accompagnying the 2 —>■ 2 parton 
scattering through the t-channel exchange are then given by 



Win) 



c 



Nc 



(13) + (24) + - {(12) + (34) + (14) + (23)} 
(13) - ^(24) + I {(12) + (34) + (14) + (23)} 



(54) 
(55) 



2Cf[{U) + (23) 

{2((12) + (34)) - (14) - (23) - (13) - (24)}], (56) 



2NcCf 
2Cf[{12) + (34) 

{2((14) + (23)) - (12) - (43) - (13) - (24)}], (57) 



2NcCf 

Considering the subprocesses i = c, . . . , g separately we obtain the ratios 

Nc 



R'/ 
R'. 



ACf 
Nc 
ACf 



(5 — cos 9 

1 , 

c 



4- 



^2vi-cose,) 



:i + cose,) + 



NcCf- 



1 

2NcCf 



(3 cos Qs — 1)- 



(58) 
(59) 
(60) 
(61) 
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At low scattering angles ^ 1 where the gluon exchange dominates all 
these ratios are 

Nr 

^1 = 7^-2.25 (^ = c,...,/) (62) 

and for the important processes c) and d) they vary only slowly with angle. 
Therefore also for the superposition ( P^) of all processes the ratio will be 
close to (|62D for not too large angles. On the other hand, these ratios differ 
markedly (by about a factor of 2) from R"f^ for the direct photoproduction. 
In the large Nc limit and at <C 1 correspond to the incoherent sum 
of two qq antennae-dipoles, while R"f^ are given by only one qq antenna. 

Therefore, when comparing the perpendicular radiation pattern in direct 
and resolved photoproduction of di-jets one expects an increased particle 
density in the latter case by about a factor of 2. It will be interesting to 
find out experimentally whether such a difference occurs already at such low 
momenta as a few hundred MeV. Such a finding would provide an important 
support for the relevance of perturbative QCD, in particular of the colour 
coherence, to the soft particle production. 



5.4 Hard processes in hadron hadron collisions 



The tests which are interesting from our point of view are very similar to 
those discussed for photoproduction in the last subsection. First there are 
the subprocesses (^)- (|4^) of parton parton collisions (at small angles) which 
we discussed in connection with resolved photoproduction. Secondly, there 
are the processes with a final state direct photon (or, equivalently a weak 
vector boson W or Z) produced together with a hadronic jet with opposite 
transverse momenta (see, e.g. 0j51|]). These processes correspond to the 
direct processes in photoproduction. The relevant subprocesses are 



a') q{pi) + q{p2) 
b') q{Pi)+9{P2) 



i{pz) + g{Pi), 
7(P3) + q{pi) 



(63) 
(64) 



and the soft bremsstrahlung is again found from 
beling 



8]) with appropriate rela- 



W (n) 



Nc 
Nc 



(14) + (24) - ^(12) 



c 



(12) + (24) - ^(14) 



N 



c 



(65) 
(66) 
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Then we obtain for the perpendicular radiation 



R 



1 



(67) 




Nc 



3 — cos©s — -^{i + COS©s) 



(68) 



c 



just as in case of direct photoproduction. 

There is an interesting difference with the case of photoproduction. At 
high primary energies the parton coUisions with initial gluons are more likely. 
Therefore in photoproduction the photon-gluon fusion process (a) without 
angular dependence in the soft particle production rate is favoured whereas 
at the hadron collider the other process {b') with angular dependence is more 
likely. This offers the possibility to study the angular dependence of the soft 
particle production in the same process. 



The analytical perturbative approach to multiparticlc production together 
with the assumption of LPHD has proven to be very successful in the de- 
scription of various inclusive characteristics of jets. It is of importance to 
investigate further the potential and limitations of this picture, in particular 
in the soft region where hadronization effects could in principle wash out the 
perturbative predictions. 

The soft particles follow a striking prediction of the perturbative anal- 
ysis. The particle density at low momentum is nearly energy independent 
over two orders of magnitude in the cms energy of e'^e" annihilation and 
this behaviour is supported now also by the data from ep collisions. More 
generally, the analytic calculations suggest a scaling law for particles at low 
transverse momentum; such a property was observed already long ago in the 
PETRA energy range. The scaling law is derived analytically for soft gluons 
which are emitted coherently from all other more energetic colour sources. 
In the perturbative description this process of single gluon emission is domi- 
nated by the Born term of 0{as) which is energy independent. At larger p± 
gluons are more frequently produced through cascading processes of higher 
order in a^. Remarkably, the hadrons follow this prediction demonstrating 
the relevance of LPHD also for soft particles. 

It is very important to further clarify whether the observed scaling be- 
haviour is just a general property of hadronization or should be considered 
as consequence of the perturbative QCD as in our picture. In the latter 
case the density of the low momentum particles should vary in a well defined 
way with the type of reactions and orientations of the primary partons in the 



6 Conclusions 
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hard process. In this paper we derive predictions for processes for which data 
are aheady available. Experimental studies of the proposed type would clar- 
ify at which scale of the low transverse momentum or energy the behaviour 
expected from the perturbative analysis sets in. 
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Figure captions 



Fig.l. a: Charged particle distribution dn/d^p as a function of particle 
momentum p. Experimental data at various cms energies p6|, p^ , ^j25|j28[] are 
compared with predictions using dn/d^p = 2-A/9- KhD^{C,E)/[4:7^E{E^ — Ql)] 
with Dl computed in MLLA from eqs. (|l],D (Qo = 0.27 GeV, Kh = 0.45), 
which approach a common limit for p —>■ 0. The detailed shape in the dashed 
region is particularly sensitive to parton level approximations, b: Same dis- 
tribution for charged pions[^,^,^,^ in comparison with predictions using 
dn/d^p = 2 ■ 4/9 ■ K'ftD^(^£;)/[47rEj|] with computed from the Limiting 
Spectrum (Qo = 0.155 GeV, Kh = 1.125). 

Fig. 2. a: The same observable as in Fig. la, but only the data collected 
at cms energies from 3 up to 58 GeV [ PB| , P^P^ , P7[| are shown, b: The same 
observable as in a, but only the data collected at LEP at cms energies a/s = 
91, 133, 161 GeV|2|,|28|j3l|,||,|3|] are shown. 

Fig. 3. a: Charged particle distribution dn/d^p at three values of momen- 
tum p = 0.2, 0.7 and 1.2 GeV as a function of cms energy ^/s. Experimen- 
tal data from e~^e~ annihilation experiment s [p^ , 2l , 32 , 27 , 2^, 2^, pl] , p^ ,^ and 
from ep experiment in the Breit frame. In the latter case two sets of data 
are shown, corresponding to different ways of sampling the data (see |^); 
in both sets data have been multiplied by a factor 2 to obtain a common 
normalization with e~^e~ data and a selection on the energy flow has been 
applied 0]. The solid line corresponds to formula (^ and the dashed one 
to formula (|lO|) . b: Distribution dn/d^p of charged kaons and of protons at 
fixed values of momentum as a function of cms energy. Experimental data 

, , and distribution at momen- 

tum p = 0.25 GeV for charged particles in e~^e~ annihilation as a function of 
cms energy ^/s. Data from |l9| , p5| , pT|p3| , p3[ . 

Fig. 5: m^dn/d^p distribution for charged pions, charged kaons, protons and 
neutral kaons in e^e~ annihilation as a function of the particle momentum 



from e^e experiments p4| , |3Tl , |30 1 . 
Fig. 4: Rapidity density at ?/ ^ ^ 



p. Experimental data from |[29|j3^j30| , |31 



Fig. 6: Theoretical prediction at ^/s = 14 GeV for the single inclusive 
particle density d'^n/dp\\d'd as a function of the angle d at four different 
values of p\\ = 0.15, 0.3, 0.5, 0.8 GeV (from right to left) according to MLLA 
and eq. ( ^2] ) (with parameters Qo = 0.5 GeV and A = 0.05). 
Fig. 7: Theoretical predictions as in Fig. 6 for d'^n/ dp\\dp^ as a function 
of transverse momentum at fixed value of = 0.8 GeV at four different 
values of cms energy ^/s = 14 (solid), 34 (long-dashed), 91 (short-dashed) 
and 172 GeV (dotted line). 
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Fig. 8: Comparison of different theoretical predictions for tfie cPn/dp\\dp±_ 
distribution as a function of transverse momentum at fixed value of p\\ — 
0.8 GeV and ^ = 34 GeV: DLA Born term (long-dashed), DLA (short- 
dashed), MLLA (solid) with Qq = 0.5 GeV and A = 0.05 and Limiting 
Spectrum (dotted line) with the same Qq and A = (normalization adjusted 
by factor — 3.5). 

Fig. 9: The spectrum of soft particles perpendicular to the event plane de- 
pends on the colour and relative orientation of the primary partons. 
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